METHOD FOR SOLVING A MULT I -GOAL PROBLEM 



The present invention relates to a method as defined in the 
preamble of claim 1. 

When the most advantageous alternative is to be selected in a 
5 situation where the final result depends on a plurality of 
factors, there often arises a conflict regarding the emphasis 
to be given to different factors. When the properties and 
ways of action of different factors are similar and commensu- 
rable, it is generally easy to develop methods in which the 
10 factors are mutually correctly weighted and the changes oc- 
curring in them are properly taken into account . 

For example, to optimize the way in which an elevator or ele- 
vator group serves a call issued by a passenger, the tradi- 
tional approach is to calculate the delays and passenger 

15 waiting times. By using coefficients, it is possible to con- 
trol the degree of importance assigned to the passenger's 
waiting time at a- floor, the passenger's traveling time in an 
elevator car and the stops during the travel of the car pro- 
posed for the passenger. As all these factors are quantities 

2 0 of time, comparing and matching them to each other will not 
involve insuperable difficulties. The goals of optimization 
can also be easily changed. 

When the factors to be optimized at the same time are not 
commensurable, it is difficult to compare them and to take 

2 5 them equally into consideration. It may be possible to accu- 

rately determine the share of individual factors in a cost 
function. However, different factors may have different de- 
grees of influence, their effects on the matter as a whole 
may appear on quite different levels, and these effects may 

3 0 even be conflicting. Thus, optimizing the cost function so as 

to reach a desired goal is a very extensive and multi- 
dimensional process. 



In the allocation of elevator calls, the objective may be to 
serve the passenger having pressed a call button as soon as 
possible and to transport the passenger to the destination 
floor without delay. On the other hand, the elevator control 

K 5 system must take into account the calls and expectations of 
other elevator passengers as well. Furthermore, the elevator 
or elevators is/are designed to take care of all internal 
transportation needs within the building, so the allocation 
of an individual call is subject to additional conditions re- 

10 lating to traffic situation, traffic intensity and available 
capacity. If the elevator control system additionally has to 
take into account the minimization of energy consumption, aim 
at reducing the number of starts of the elevator or park any 
elevators that may be free in the current traffic situation 

15 at certain floors by considering overall advantages, then 
managing the cost function by prior-art methods is an impos- 
sible task. 

The object of the invention is to disclose a new method for 
optimizing a solution to a problem situation in which the so- 
20 lution is influenced by a plurality of factors that are not 
commensurable quantities. To achieve this, the method of the 
invention is characterized by the features presented in the 
characterization part of claim. 

By the method of the invention, a multi-goal optimization 
25 problem can be solved quickly and reliably so that different 
factors contributing to the optimization are weighted in a 
desired manner. The computation time needed in the optimiza- 
tion can be limited to a short time so that, in situations 
where the computing time is limited, alternative solutions 
30 are considered when a decision is being made. E.g. in eleva- 
tor group control applications, in which allocation decisions 
have to be made repeatedly and for constantly changing cost 
functions, speed and efficiency are of primary importance. 



By utilizing the properties of genetic algorithms, sub- 
functions and overall optimization can be executed advanta- 
geously and very quickly with reasonable computing capacity. 

r In the following, the invention will be described in detail 
5 by the aid of an example of its embodiments with reference to 
the attached drawings, wherein 

- Fig. 1 visualizes a multi-goal optimization problem 

- Fig. 2 represents the differences between the distribu- 
tions of the goals of the multi-goal problem 

10 - Fig. 3 illustrates an approach according to the invention 

- Fig. 4 represents normalized distributions of cost func- 
tions 

- Fig. 5 presents an example based on a genetic algorithm 
according to the invention. 

15 In the following, a solution to a multi-goal problem is de- 
scribed where the objectives are, on the one hand, optimiza- 
tion of energy consumption and, on the other hand, optimiza- 
tion of passengers' call times. In mathematical terms, the 
optimization problem for solution alternative A of the total 

20 cost function J can be expressed by the equation 

J (A) = EWiCMA) , 

Where Ci represents an individual cost function, in 
this example call time and energy consumption for 
alternative A and 



25 



Wi represents a weighting coefficient assigned to 
the individual cost function. 



In this case, the solution to the optimization problem is 
minimization of function J. A problematic question is how to 
define correct values for the weighting coefficients. If a 
given cost function, such as call time, gets a high weight- 
5 ing, then it will become dominating and the influence of the 
other factors will remain marginal. Also, a small cost func- 
tion may have a very small influence. 

Referring to Fig. 1 and 2, let us consider the optimization 
of passengers' call times and energy consumption of the ele- 
vator in the same space A c of allocation solutions (reference 
number 1) , which contains all possible solutions for serving 
the calls active in the elevator group. The allocation alter- 
natives can be divided into two sub- spaces CT (2) and E (3) 
according to their relation to call times on the one hand and 
to energy consumption on the other hand. These spaces have 
statistical properties such as distribution, expectation 
value £ and variance a 2 . The statistical properties of these 
two spaces are described in Fig. 2. In addition to the dif- 
ference of units of measurement - the unit for call time is 
second while the unit for energy consumption is Joule - the 
quantities also differ from each other in respect of statis- 
tical properties, as appears from Fig. 2. 

Besides being non- commensurable, the targets of optimization 
are also to be weighted in different ways in different situa- 

2 5 tions. For example, the task may be to find a solution in 

which energy consumption has a weight of 30 % and call times 
have a weight of 70 %. 

Theoretically, normalized cost factors % can be defined if 
the expectation value £ and variance a 2 of the cost space are 

3 0 known, by the equation 

X = (C - I) /a. 



10 



15 



20 



5 



In practical solutions, such a procedure is not viable be- 
cause going through the entire space to be considered is a 
, task too laborious and in most cases impossible. Instead, the 
5 expectation value and variance can be approximated by using 
their sample equivalents, sample average \i and sample vari- 
ance s 2 . The normalized cost function can thus be expressed 
in the form 

X « (C - \i) /s. 

10 The sample average |i is normally distributed with variance 
<T /n, which can well be used to estimate the required number 
of samples n. Fig. 3 presents a drawing visualizing the 
utilization of a sample in the definition of normalized func- 
tions. Where applicable, the designations and reference num- 

15 bers used in Fig. 3 are the same as in Fig. 1. From sub-space 
2, a sample 12 has been taken, which contains a certain set 
of the elements of space 2. In the example of allocation of 
elevator calls implemented using a genetic algorithm that is 
described below, this set of samples preferably consists of 

20 members of a first generation of solutions. In a correspond- 
ing manner, a sample 13 has been taken of sub-space 3. For 
the samples depicted in Fig. 3, the statistical quantities 
sample average fi and sample variance s 2 are defined, which 
approximately describe the statistical quantities expectation 

25 value £ and variance a 2 for the entire sub-spaces 2 and 3 in 
the manner described above. 

Fig. 4 visualizes the relationship between the normalized 
cost functions. As the cost functions are commensurable, they 
can be added together and their sums can be evaluated by the 
30 same criteria. As indicated in Fig. 4, the normalized cost 
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function obtained for call time is CT = (CT - |i C t) /s CT and 
correspondingly the normalized cost function for energy con- 
sumption is E = (E - |LIe)/Se. The normalized total cost func- 
tion, which is to be minimized, is correspondingly 

5 J = KctCT + K E E, 

where Kc T and K E are drive-specific coefficients to 
be determined separately. 

In the following embodiment example, the implementation of 
multi-goal optimization using a genetic algorithm is de- 
10 scribed. Below is a short summary of the application of a ge- 
netic algorithm to the allocation of elevator calls. For a 
more detailed description, reference is made e.g. to patent 
specification US 5932852. 

When calls are allocated by means of a genetic algorithm, 
15 each landing call is encoded as a gene of a call chromosome. 
The position of the gene in the chromosome represents an ac- 
tive landing call, and correspondingly the value of the gene 
represents the elevator car proposed to serve the landing 
call. Each chromosome represents one alternative solution to 

2 0 the allocation problem that is able to serve the active 

calls. From the chromosomes, a population typically compris- 
ing about 50 chromosomes or solution alternatives is formed. 
For each chromosome in the population is determined a so- 
called Fitness value, which consists of the sum of the cost 
25 functions of the elevators serving active calls. The cost 
functions are defined on the basis of selected criteria, and 
their values are computed using a model of each elevator. 

After the Fitness values of all the chromosomes have been de- 
termined, they are listed in order of Fitness values. From 

3 0 the chromosomes, new generations are formed by genetic algo- 
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rithm methods. After about 20-50 generations, the best al- 
ternative can be found, and this alternative is selected to 
serve the active landing calls. 

/ Fig. 5 visualizes an example embodiment of the invention in 
5 which a multi-goal problem is solved by utilizing both nor- 
malization of non-commensurable cost functions and methods of 
allocation based on a genetic algorithm. As for the formation 
of chromosomes and computation of the Fitness values, refer- 
ence is made to patent specification US 5932852. 

10 On the basis of the active landing calls and car calls, the 
chromosomes 40 of the first population are generated, on the 
basis of which the Fitness values of the allocation alterna- 
tives corresponding to the chromosomes are determined, con- 
sidering both call time optimization CT and energy consump- 

15 tion E, in a computation unit 42. In the example presented in 
Fig. 5, the elevator group comprises two elevators, elevator 
A and elevator B. For each elevator, an elevator model 44 and 
46, respectively, has been formed, these models comprising 
the required elevator-specific information for the calcula- 

20 tion of the cost functions. Based on this information and the 
active calls to be served, cost functions are determined in 
the computation unit for both call times CT A and CT B and en- 
ergy consumption E A and E B . A cost function CT for the call 
times of the entire elevator group for a given allocation al- 

25 ternative is obtained as the sum CT = CT A + CT B , and a cost 
function E for energy consumption in the entire elevator 
group is obtained correspondingly from the sum E = E A + E B . 
These partial cost functions for call times and energy con- 
sumption are stored in tables 48 and 50 of partial Fitness 

3 0 values. 

A first population is produced e.g. in the manner described 
in patent specification US 5932852. Based on the partial Fit- 



ness values of this first population, i.e. on the values of 
the partial cost functions, sample averages |A PF i and jj. pf2 and 
sample variances s 2 PF i and s 2 PF 2 for a sample according to the 
first population are determined in the manner specified in 

- 5 Fig. 3 and formulas 1-3. These sample quantities jx and s 2 are 
used in the calculation of the Fitness value 54 of a chromo- 
some. In the determination of the Fitness value, a weighting 
coefficient K PF1 and K PF2 (block 58) defined for the partial 
cost function by the operator 56, e.g. the owner of the 

10 building, is taken into account. The calculated results con- 
stitute the total Fitness value of the chromosome and they 
are stored in a table 60. On the basis of these values, the 
best solution alternatives of the population are evaluated. 
In the next populations, the sample quantities \x and s 2 are 

15 utilized, which are used to normalize the partial cost func- 
tions, whereas the other factors used a basis of calculation 
change in a manner determined by the genes of the chromosome 
and the elevator models. 

In the embodiment example presented in Fig. 5, the normaliza- 
2 0 tion of the partial cost functions and the calculation of the 
values of the normalized cost functions are performed in 
block 54, whereas the calculation of the values of the sub- 
functions, in this case call times and energy consumption, is 
performed in block 45, taking the call situations and eleva- 
25 tor models into account. 
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